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We present a spin-1/2 bilayer model for the quantum 
order-disorder transition which (i) can be solved by mean-field 
theory for bulk quantities, (ii) becomes critical at the tran- 
sition, and (iii) allows to include intralayer frustration. We 
present numerical data (for systems with up to 240 sites) and 
analytical results for the critical coupling strength, ground- 
state energy, order parameter and for the gap. We show that 
the critical coupling decreases linearly with frustration. 
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The model of two coupled antiferromagnetic layers, the 
'bilayer model', has attracted much attention in the last 
years [jl] || . The physics of this model is dominated by 
the competition between the in-plane coupling Ji | and the 
inter-plane coupling J\2- A transition occurs from an or- 
dered antiferromagnetic state to a spin-liquid state with 
a gap at a certain critical ratio J c of Jnj J\\- Millis and 
Monien Q proposed that certain anomalies in the mag- 
netic response of Y BaiCu-sO^+x, interpreted commonly 
as evidence for a spin gap might be explained within the 
bilayer model. 

Solvable model Hamiltonians play an important role 
in condensed matter theory, as they provide new insights 
and because they can be used as references for approxi- 
mate theories. Here we propose a non-trivial 'long-range' 
version of the bilayer model for which we are able to solve 
for all quantities of interest, while retaining the crucial 
ingredients driving the quantum order-disorder transi- 
tion. In this model the in-plane coupling Ji is of long- 
range Lieb-Mattis type while the inter-plane coupling, 
J12, remains short ranged. As a result, bulk expectation 
values for this system are accessible by analytical con- 
siderations, we are able to determine the exact ground- 
state wavfunction up to intensive corrections. Excita- 
tions, remain non-trivial however, but the simplicity of 
the in-plane interaction facilitates the numerical analy- 
sis to obtain excellent control over finite size corrections. 
We determine the functional dependence of the gap as 
a function of J12/J1, presenting the first rigorous study 
of the impact of in-plane frustration on the stability of 
the spin-liquid phase. We find that the critical J c is re- 
duced linearly with the in-plane frustration |^||. This 
result is of particular interest since a substantial inter- 
layer coupling constant J12 was found in recent infrared 



absorption experiments on YiBa2Cu306 R and by a de- 
tailed quantum chemical calculation ||||]. 

Model: We define the 'long-range' bilayer Hamiltonian 
as H LR = Y J ~« H i+ #12, with 



H 1 = (Jx/N) S A7 • S B7 
H12 = J12 JJi • Sj^, 
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and J\, J2, J12 > 0. Here S;. 7 are the spin-1/2 operators 
of the first (7=1) and of the second (7 = 2) layer 
respectively and S A -y = X lG A S ^7' S b 7 = X^gb S *,7 
are the total-spin operators of the A/B sublattice in the 
respective layers. Each sublattice contains N spins, the 
total number of sites is N s — AN. 

The first term of the intra-layer interaction (Ji), 
favours an ant i- alignment of the sublattice spins, i.e. an 
antiferromagnetically ordered state. The second, frus- 
trating, intra-layer interaction (J2) favours minimal sub- 
lattice magnetisation and therefore tends to suspress the 
antiferromagnetic order. The local inter-plane coupling 
term H12 favours the formation of local singlets. As a 
result, this model contains the relevant interactions to 
exhibit a quantum order-disorder transition as the rel- 
ative strength of intra- and inter-layer interactions are 
varied. For small values of the interlayer coupling J12 
the ground-state is characterized by long-range antifer- 
romagnetic order and gapless magnetic excitations. The 
spin- liquid state, realized for large J12, is characterized 
by a finite spin gap and the absence of antiferromagnetic 
order. 

The long-range bilayer Hamiltonian Hlr has a close 
relation to its 'short-range' counter-part, Hsr = 
^HS R + H 12 , with 



(2) 



<hj> 



where the symbols < i.j > and [i, I] denote pairs of n.n. 
and n.n.n sites on a square lattice respectively and where 
the terms ~ J|| and ~ JL correspond to the term ~ J\ 
and ~ J2 in ([!]) respectively. In order to make connection 
between Ulr and H$r we compare the energies of the 
Neel state in both models for the case J 2 = J'u = 0. This 
introduces a correspondence relation Ji=4J| between 
the respective couplings. In Hsr, the order-disorder tran- 
sition occurs at J\2/J\\ ~ 2.51 0. We expect the true 
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critical interlayer coupling for Hlr to be somewhat larger 
than the value J12/J1 ~ 2.51/4 sa 0.628 suggested by this 
naive rescaling, since intra-plane quantum fluctuations 
are suspressed by the long-range nature of J\. 

Results: In the limit of decoupled layers J12 = the 
expectation value of the energy per layer is 
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2J 2 - J 



2N 



-[S A1 (S A1 + 1) + S B1 (S B1 + 1)}. 



where Sai> Sbi are the eigenvalues of Sai, S^i and of the 
total-spin operator Sai + S^i respectively [ ^0|JlT| . For 
J 2 < Ji/2 the ground state is a singlet with 5i = and 
maximal sublattice magnetization, Sai = Sbi = N/2. 
In the thermodynamic limit the ground-state energy per 
site is then (2J 2 — Ji)/8. In the opposite limit J 12 = 00 
the ground state is a product of interlayer singlet pairs 
with energy — 3/8J 2 per site and a gap of order J 12 to 
the lowest excited state, a spin triplet. 

Variational Analysis: In the presence of a finite inter- 
layer coupling a natural ansatz for a variational wave- 
function interpolates between the two extremal regimes: 

> = - a ^.A + S^Sr 2 ]) \AF >, (3) 



where \AF > is the Neel state |l!q ]. A straightforward 
minimization of the energy yields the critical coupling 
constant J c = J\ — 2 J 2 . The optimal value for the varia- 
tional parameter is found to be a = 1, for J 12 > J c , and 
to be a = J c /Ji 2 - \J{J c lJ\i) 2 - 1 for J 12 < J c - 

Of interest are the expectation values of the ground- 
state energy per site: 



e N =< H > /(AN) + 3/8J12, 



(4) 



where we have subtracted the energy of the product state 
realized for J 12 — > 00, and for the in-plane staggered 
order parameter, 



on 



V\<S A1 -S B1 + S A2 -S B2 >\/{AN). (5) 

The variational ground-state energy per site (^) and the 
order parameter (^|) are then computed to be 



e = -^[l- J12/JJ 2 



0(1 



VI - (Jl2/Jcf 



(6) 
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for J 12 < Jc- For Ji 2 > J c both eo and 00 are zero. For 
■^12 > Jc the ground state is a product of inter-layer spin 
singlets, for J\ 2 < Jc an ordered antiferromagnet. 

The quality of a trial wavefunction I'J' > is measured 
by its variance, 



< $\(H- < H >) 2 |* > 



(8) 



which is zero for an exact eigenstate of H . In general, 
variational wavefunctions which are not exact cigcnstates 



lead to an extensive variance, a ~ O(N) JfJ[. For the 
wavefunction given by (^), however, we find an intensive 
variance, cr(a) ~ O(N ), which indicates that the ex- 
pressions for the variational ground-state energy eo, Eq. 
^, and order parameter, Eq. ^, are exact in the thermo- 
dynamic limit. However, since the energy is known only 
up intensive correction, the gap cannot be determined in 
this framework. 

Numerical Analysis: The Hamiltonian (0) is separately 
invariant under any permutation within the A-sublattice 
and within the B-sublattice. The ground-state singlet 
and the lowest excited triplet state belong to the symmet- 
ric subspace of the permutation symmetry, as we found 
by diagonalizing small clusters with N s =4 — 20 spins 
within the full Hilbert space. We conjecture that the 
symmetry of the lowest singlet state and of the lowest 
triplet state does not change for N s > 20. Within the 
symmetric representation the size of the Hilbert space 
grows only ~ N^, i.e. algebraically. We have numerically 
diagonalized systems with up to N s = AN = 240 sites, 
using a specially adapted exact diagonalization technique 
Jl4j]. The results for the ground-state energy per site, are 
presented in Fig. for J2/J1 = 0, 0.2. The width of 
the transition from the ordered state at small J\ 2 / J\ to 
the quantum disordered spin liquid at large Jyi/ J\ is at- 
tributable to finite size effects US 

Using finite-size scaling we have estimated the critical 
exponents of the finite-size corrections to the data for 
the order parameter presented in Fig. [|. We found the 
finite-size corrections to scale like TV" 1 , iV" 1 / 3 and TV -1 / 2 
for J12 < Ji, J12 = Ji and Jyi > Ji (with J 2 = 0) 
respectively. A similar behaviour has been found for the 
long range Ising model in a transversal field |l6| ]. 

We have extrapolated the numerical data obtained 
from clusters with up to 240 sites to the thermodynamic 
limit via the formula limjv^oo ejv = e x + B /N + C/N 2 + 
. . ., where we have used the 1/N finite size scaling of the 
Lieb-Mattis antiferromagnet as the leading term |]Io| , |rT| . 
The 1/N scaling for the finite-size corrections of the 
ground-state energy is also found in a least square fit to 
the numerical data presented in Fig. [I]. For all parameter 
values considered, i.e for J12/J1 — 0, 0.1, . ..,1.0 and for 
J 2 /Ji = 0, 0.125, . . .0.5, we found agreement to at least 
six digits between the such obtained eoo and the analyti- 
cal formula for e a (||) . This very good agreement between 
the numerical and analytical results gives us confidence 
that the gap, which is accessible only to numerical calcu- 
lation can be determined with similar precision. 

In Fig. [2] we present the numerical results of the gap, 
Ajv, for clusters with up to N s — 240 sites and J 2 = 
together with the mean-field prediction for the gap, 



A 



Jc + Jl2 



MF = 



6{J C -Jii) + J12 0(J 12 - Jc)- (9) 



Due to the spin-rotational invariance of (0) the gap Ajv 
differs non-trivially from the mean-field result. Note, 
that this is possible only because the variance (||) of the 
mean-field solution (||) is intensive and not zero. Fig. || 
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also contains a conjecture for the gap in the thermody- 
namic limit. 



Jl2\/1-Jc/Jx 



(10) 

for J12 > J c = Ji — 2J2. In the insert of Fig. ^ we have 
plotted log(AAr — A Q ) versus log(l/7V). We note that the 
data fall on a straight line with slope one for all J12 > 
J c , suggesting that ( jToj ) is indeed the correct expression 
for the gap in the thermodynamic limit. At the critical 
point, J12 = J 1 — 2 J 2 , the finite-size corrections scale like 
{1/N) V . We estimate numerically v — 0.3, a value which 
deviates from mean-field predictions and for which we do 
not have an explanation at present. 

Discussion: In (|^) we defined the order parameter on 
of the bilayer system. It is of interest to examine with 
5 N = yj\ < S A1 ■ S B2 + S B i ■ Sa2 > I / (42V) the rel- 
ative ordering of the two layers. An inspection of the 
variational wavefunction (^|) shows that on and on are 
directly proportional in the thermodynamic limit. A sim- 
ilar result has been found in a numerical study of the 
short range model [Q. The equivalence of on and Ojv 
implies that the magnetic order parameters of the two 
layers are aligned for any < J12 < J c - This type of 
magnetic correlation has been found in neutron scatter- 
ing experiments on YBaCuO system jlTj . Physically we 
understand that any J12 > lifts the degeneracy of the 
two intra-layer Neel states and correlates the two layers 
magnetically. 

In the long range bilayer model the order-disorder 
transition is second order with mean-field exponents for 
ground-state properties |l8). The gap opens continu- 
ously on the disordered side, with an exponent for the 
finite-size correction which deviates from the predictions 
of mean-field theory. While the in-plane quantum fluctu- 
ations are too weak to influence bulk properties like the 
ground-state energy, they arc important for the gap be- 
tween the exact ground-state, which is an overall spin sin- 
glet, and the lowest excited triplet state. In the asymp- 
totic limit J12/J1 — > 00 the gap ( |lC)|) is J12 but for any 
J12 1 J\ < 00 the gap is reduced by the in- plane quantum 
fluctuations and is strictly smaller than J12. 

Our rigorous result that the critical coupling strength 
depends linearly on the frustration is consistent with a re- 
cent effective-action mean-field approach |19 to the frus- 
trated bilayer model 0. Furthermore this results does 
support, to some extend, the possible explanation of the 
observed spin-gap in the YBaCuO System, in view of 
the recent estimate of J12/J1 ~ 0.55 by infrared absorp- 
tion Q, by a doping- frustrated bilayer system Let 
us also note in this context that the linear reduction of 
the sublattice magnetization by the inter-layer coupling 
(0) is due to quantum fluctuations only and that the 
order-disorder transition found in Hlr has no classical 
analagon. Eq. (Q) represents therefore the first rigorous 
result for the reduction of a classical order paramter by 
quantum fluctuations in dimensions greater than one. 

In summary we have presented a novel, solvable model 
for the quantum order-disorder transition. The quantity 



of interest, the spin gap, is critical at the transition while 
the bulk properties of this model can be determined an- 
alytically. The intra-layer quantum fluctuations are only 
intensive, which allows a complete characterization of the 
qunatum order-disorder transition. The inter-layer quan- 
tum fluctuations remain extensive, introducing the criti- 
cality neccessary for the transition. 
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FIG. 1. For N = 10,20,30,40,50,60 and a) J 2 /Ji = 
and b) J2JJ1 = 0.2 the numerical results for the order pa- 
rameter (H) (open squares, smallest N at the top) and for the 
ground-state energy per site (^) in units of J\ (open circles, 
smallest N at the bottom). The dashed lines are the analytic 
expressions corresponding to Eq. (0) and Eq. (M). 



FIG. 2. For N = 10,20,30,40,50,60 and J 2 /J\ = 
the numerical results for the gap Ajv (open circles, small- 
est N at the top). The dashed line is the analytic ex- 
pression corresponding to Eq. (|lo|), the dotted line is the 
mean-field prediction, Eq. (|^). Insert: A log-log of the data 
for J12/J1 = 1.0, 1.2, 1.6, 2.0. 
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